import math
from matplotlib import pyplot

gamma0=0.025
sigma0=[1.19155,1.52364,2.1085,3.07389,3.85139,4.21462,4.36921,4.59828,4.56225,4.66273,4.82263,4.87958]
delta=[1.18813,1.22498,1.1991,1.14654,1.01877,0.880983,0.739985,0.606701,0.482742,0.285949,0.104621,4.171e-12]
omega=[0.5,1,2,5,10,15,20,25,30,40,50,100]

f = open("G1G2Gamm0.025FREQ.dat","w")
for i in range(len(sigma0)):
    G1=sigma0[i]*math.cos(delta[i])/gamma0
    G2=sigma0[i]*math.sin(delta[i])/gamma0
    f.write(str(omega[i])+' '+str(G1)+' '+str(G2)+' '+str(sigma0[i])+' '+str(delta[i])+'\n')
    print "G': ",G1, "G'': ",G2
f.close()

